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ABSTRACT 

We examine the effect of weak, small scale magnetic field structure on the 
rate of reconnection in a strongly magnetized plasma. This affects the rate of 
reconnection by reducing the transverse scale for reconnection flows, and by 
allowing many independent flux reconnection events to occur simultaneously. 
Allowing only for the first effect and using Goldreich and Sridhar's model of 
strong turbulence in a magnetized plasma with negligible intermittency, we 
find a lower limit for the reconnection speed ~ VaR^^ U '.M 3//4 , where Va is the 
Alfven speed, TZl is the Lundquist number, and M. is the large scale magnetic 
Mach number of the turbulence. We derive an upper limit of ~ VaM 2 by 
invoking both effects. We argue that generic reconnection in turbulent plasmas 
will normally occur at close to this upper limit. The fraction of magnetic energy 
that goes directly into electron heating scales as IZ^ 

2/5^8/5 

and the thickness 

of the current sheet scales as TZ L 3 ^Ai~ 2 ^ 5 . A significant fraction of the magnetic 
energy goes into high frequency Alfven waves. The angle between adjacent field 
lines on the same side of the reconnection layer is ~ TZJ^^ 5 A4 6 ^ 5 on the scale 
of the current sheet thickness. We claim that the qualitative sense of these 
conclusions, that reconnection is fast even though current sheets are narrow, is 
almost independent of the local physics of reconnection and the nature of the 
turbulent cascade. As the consequence of this the Galactic and Solar dynamos 
are generically fast, i.e. do not depend on the plasma resistivity. 

Subject headings: Magnetic fields; Galaxies: magnetic fields, ISM: molecular 
clouds, magnetic fields 
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1. Introduction 

Magnetic reconnection is a long standing problem in theoretical magnetohydrodynamics 
(MHD). This problem is closely related to the hotly debated problem of the magnetic 
dynamo (see Parker 1979, Moffatt 1978, Krause & Radler 1980). Indeed, it is impossible 
to understand the origin and evolution of large scale magnetic fields without a knowledge 
of the mobility and reconnection of magnetic fields. In a typical astrophysical plasma, 
resistivity is very small and flux freezing should be an excellent approximation to the 
motion of magnetic field, while dynamo action invokes a constantly changing magnetic field 
topology [] This means an astrophysical dynamo requires efficient reconnection despite very 
slow Ohmic diffusion rates. 

The literature on magnetic reconnection is rich and vast (see, for example, Biskamp 
1993 and references therein). In the simplest model of magnetic reconnection ( Parkerj 
1957| , (Sweet 1958| , see §3.1) the reconnection speed lies below the Alfven speed by a factor 



~ Tl^ 2 = (ti/VaL x ) 1 I 2 , where rj is the resistivity, Va is the Alfven speed, and L x is the 
length of the current sheet, assumed to be determined by the large scale geometry of the 
problem. The large ratio of L x to the width of the current sheet, A, equal to the ratio of 
the resistivity to the reconnection speed, drives the reconnection speed towards zero for 
highly conducting fluids. In particular, this is many orders of magnitude too small for 
astrophysical reconnection. Consequently, there has been intense interest in models which 
result in higher reconnection speeds. 

In general, we can divide schemes for fast reconnection into those which alter the 
microscopic resistivity, broadening the current sheet, and those which change the global 
geometry, thereby reducing L x . An example of the latter is the suggestion by Petschek 
(1964) that reconnecting magnetic fields would tend to form structures whose typical size 
in all directions is determined by the resistivity ('X-point' reconnection). This results 
in a reconnection speed of order V^/ln"^. On the other hand, attempts to produce 



such structures in simulations of reconnection have been disappointing QBiskamp 1984] , 



Biskamp 1986|) . In such simulations the X-point region typically collapses towards the 



Sweet-Parker geometry as the Lundquist number becomes large (Biskamp 1996| , Wang, Ma 



& Bhattacharjee 1996, Wang & Bhattacharjee 1996). One way to understand this collapse 
is to consider perturbations of the original X-point geometry. In order to maintain this 
geometry reconnection has to be fast, which requires shocks in the original (Petschek) 



1 Merely winding up a magnetic field can increase the magnetic field energy, but cannot increase the 
magnetic field flux. We understand the dynamo in the latter sense. The Zel'dovich "fast" dynamo (Vainshtcin 
& Zel'dovich 1972) also invokes reconnection for continuous dynamo action (Vainshtein 1970). 
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version of this model. These shocks are, in turn, supported by the flows driven by fast 
reconnection, and fade if L x increases. Naturally, the dynamical range for which the 
existence of such shocks is possible depends on the Lundquist number and shrinks when 
fluid conductivity increases. The apparent conclusion is that at least in the collisional 
regime reconnection occurs through narrow current sheets. | 

In the collisionless regime the width of the current sheets may be determined by the ion 
cyclotron (or Larmor) radius ( Parker 1979|) or by the ion skin depth ( [Ma fc Bhattacharjee 



|1996| , |Biskamp, Schwarz, fc Drake 1997| , [Shay, Drake, Denton, fc Biskamp 1998|) . This 



often leads to a current sheet thickness which is much larger than expected ('anomalous 
resistivity'). Furthermore, anomalous resistivity can be combined with X-point reconnection 
with the result that the X-point is somewhat less unstable to collapse. Nevertheless, while 
this effect can dramatically increase reconnection speeds in simulations and in laboratory 
experiments it does not change the qualitative conclusion that narrow current sheets give 
reconnection speeds in astrophysical plasmas that are unacceptably low. 

In order to explain why astrophysical magnetic fields do not reverse on very small 
scales, researchers have usually appealed to an ad hoc diffusivity which is many orders 
of magnitude greater than Ohmic diffusivity. This diffusivity is assumed to be roughly 
equal to the local turbulent diffusion coefficient. While superficially reasonable, this choice 
implies that a dynamically significant magnetic field diffuses through a highly conducting 
plasma in much the same way as a tracer of no dynamical importance 0. This is referred 
to as turbulent diffusivity and denoted rjt, as opposed to the Ohmic diffusivity r\. Its name 
suggests that turbulent motions subject the field to kinematic swirling and mixing. As 
the field becomes intermittent and intermixed it can be assumed to undergo dissipation at 
arbitrarily high speeds. 

Parker (1992) showed convincingly that the concept of turbulent diffusion is ill-founded. 
He pointed out that turbulent motions are strongly constrained by magnetic tension and 
large scale magnetic fields prevent hydrodynamic motions from mixing magnetic field 



2 A very different approach explored by Vishniac (1995a, b) appealed to the concentration of magnetic field 
in narrow flux tubes within a high (3 plasma, thereby reducing L x to the flux tube radius. Given the small 
filling factor of the magnetic field and the high Alfven velocity inside the flux tubes it is possible to obtain 
fast reconnection even using the Sweet-Parker geometry. However, flux tube formation requires initially high 
reconnection rates. 

3 Reconnection in the presence of ambipolar diffusion discussed in Vishniac & Lazarian (1998) vaguely fits 
into the category of "enhanced A" . Indeed, the neutrals do not feel magnetic field and leave the reconnection 
zone forming an outflow that is limited by ambipolar diffusion rather than Ohmic diffusivity. 

4 In the case of galaxies the ratio of the two diffusivities is ~ 10 20 . 
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regions of opposing polarity unless they are precisely anti-parallel. Both numerical and 
analytic studies (see |Cattaneo fc Vainshtein 1991| , [Kulsrud fc Anderson 199"2| , |Gruzinov fc 
piamond 1994[ ) confirm that the traditional flRuzmaikin et al. 1988| ) theory of kinematic 
dynamos is seriously and fundamentally flawed. On the other hand, observations of the 
solar corona and chromosphere seem to show that reconnection often takes place at speeds 
of ~ O.IVa (cf. Dere 1996, Innes et al. 1997 and references contained therein). Evidently at 
least some astrophysical plasmas can undergo reconnection on short time scales, regardless 
of the apparent theoretical difficulties involved. 

An associated problem for dynamo theories is the growth of small scale random 
magnetic fields. Assuming magnetic fields are allowed to intermix freely, as the turbulent 
diffusion model requires, Low (1972) and Krause & Radler (1980) estimated that the ratio 
of the mean field to the random field is expected to be of the order \Jrj t /ri 3> 1. At the same 
time, observations indicate that these field components are not very different in galaxies 
at the present epoch (Beck et al. 1996). Furthermore, Kulsrud and Anderson (1992) have 
suggested that the growth of a random field should suppress dynamo action. 

This unsatisfactory state of affairs has stimulated recent attempts to find new 
approaches to the problem of the astrophysical dynamcQ (e.g. Parker 1992, Parker 1993, 
Vishniac 1995a,b, Kulsrud et al 1997) that have deficiencies which we will discuss in more 
detail below. 

Below we explore a different approach to the problem of rapid reconnection. This 
paper deals with the issue of magnetic reconnection^ in the presence of a weak random 
field component. Magnetohydrodynamic turbulence guarantees the presence of such a field, 
although its amplitude and structure clearly depends on the model we adopt for MHD 
turbulence, as well as the specific environment of the field. It is also obvious from the above 
discussion that as long as the kinematic approximation holds, that is as long as magnetic 
forces are negligible, Q the random component of the magnetic field will quickly come to 
dominate the large scale one. We will assume throughout that we can describe resistive 



5 The astrophysical dynamo is just an example of a process whose dynamics is determined by reconnection 
rates. Solar flaring (see Dere 1996) and star formation in molecular clouds (for example see Lubow & Pringle 
1996) are other examples. 

6 The mode of reconnection discussed here is sometimes is called free reconnection as opposed to forced 
reconnection. Wang & Bhattacharjee (1992) define free reconnection as a process caused by a nonideal 
instability driven by the free energy stored in an equilibrium. If the equilibrium is stable, reconnection can 
be forced if a perturbation is applied externally. 

7 This approximation should be valid, at least, at the initial stages of magnetic field generation. 
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effects with an Ohmic resistivity rj. However, our qualitative conclusions are not sensitive 
to this choice. We will show that if we allow for the more complicated geometry expected 
from a turbulent cascade, the speed of reconnection is fast, independent of the plasma 
resistivity, and a function only of the overall level of turbulence. Although we will adopt a 
specific model of MHD turbulence in this work, we will also show that our conclusions are 
not sensitive to this choice. 

There are two mechanisms which lead to this result. First, if we consider two adjoining 
regions with very different large scale fields, the small scale structure of both fields implies 
that the typical scale for a reconnection event is much smaller than the overall size of the 
system, and goes to zero in the limit of an ideal fluid. Second, since field lines wander in 
and out of the reconnection zone, the overall rate of flux reconnection can be much greater 
than the rate for a single event. We note that it is usually believed that for reconnection to 
be rapid in the limit of rj — > a current singularity should develop (Park et al 1984). Our 
model does not require such singularities. Indeed, we show that while the amount of Ohmic 
dissipation tends to as r\ — > 0, the magnetic field has a weak stochastic substructure on 
the smallest scales and the rate of the flux reconnection does not decrease. 

In §2 of this paper we will introduce the scheme of reconnection in the presence of 
stochastic field, while in §3 of this paper we will describe the results of a simple model 
for the small scale structure of a large scale magnetic field embedded in a turbulent 
medium. A fuller discussion of this model is given in Appendix A. In §4 we will estimate 
the reconnection speed in the presence of small amounts of stochasticity, and apply these 
estimates to the model described in §3. In §5 we describe some of the implications of this 
work. We give our conclusions in §6. 

2. The Problem 

Sweet-Parker reconnection (Sweet 1958, Parker 1957) provides a minimum value for 
reconnection rates. In this model oppositely directed magnetic fields are brought into 
contact over a region of size L x . Magnetic fields reconnect along a very thin Ohmic diffusion 
layer L y ss r)/V rec and fluid is ejected from this layer at a velocity of the order Va in a 
direction parallel to the local field lines. The layer in which Ohmic diffusion takes place is 
usually referred to as the current sheet. Here we will refer to the volume where the mean 
magnetic field strength drops significantly as the reconnection zone, in order to allow for 
the presence of collective effects which may broaden the reconnection zone well beyond 
the current sheet. The reconnection velocity in the Sweet-Parker picture is determined by 
the constraint imposed by the conservation of mass condition V rec L x VaL v . Although 
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this model is two dimensional, it can be generalized to three dimensions by allowing the 
two magnetic field regions to share a common field component, which has the effect of 
rotating them so that they are no longer exactly anti-parallel. This has no effect on the 
Sweet-Parker reconnection process (see Fig. 1). However, it does change the nature of the 
constraint somewhat. In addition to ejecting matter from the reconnection zone, we must 
also allow for the ejection of the magnetic flux due to the common field component. This 
is, in effect, the same constraint in this case. 

We consider the case in which there exists a large scale, well-ordered magnetic field, 
of the kind that is normally used as a starting point for discussions of reconnection. This 
field may, or may not, be ordered on the largest conceivable scales. However, we will 
consider scales smaller than the typical radius of curvature of the magnetic field lines, or 
alternatively, scales below the peak in the power spectrum of the magnetic field, so that 
the direction of the unperturbed magnetic field is a reasonably well defined concept. In 
addition, we expect that the field has some small scale 'wandering' of the field lines. On 
any given scale the typical angle by which field lines differ from their neighbors is <C 1, 
and this angle persists for a distance along the field lines Ay with a correlation distance A^ 
across field lines. 

The modification of the mass conservation constraint in the presence of stochastic 
magnetic field component is self-evident. Instead of being squeezed from a layer whose 
width is determined by Ohmic diffusion, the plasma may diffuse through a much broader 
layer, L y ~ (y 2 ) 1 ^ 2 (see Fig. 2), determined by the diffusion of magnetic field lines. The 
value of (y 2 ) 1 / 2 can be determined once a particular model of turbulence is adopted (see §3), 
but it is obvious from the very beginning that this value is determined by field wandering 
rather than Ohmic diffusion as in the Sweet-Parker case. 

In the presence of a stochastic field component, magnetic reconnection dissipates field 
lines not over their entire length ~ L x but only over a scale Am <C L x (see Fig. 2b), which is 
the scale over which magnetic field line deviates from its original direction by the thickness 
of the Ohmic diffusion layer Aj 1 « v/^rec,iocai- If the angle of field deviation does not 
depend on the scale, the local reconnection velocity would be ~ Va4> and would not depend 
on resistivity. We claim in §3 that does depend on scale. Therefore the local reconnection 
rate Vrec 

local is given by the usual Sweet-Parker formulae but with An instead of L x , i.e. 
Vrec,iocai ~ Va(V/iA||/ 'v) -1 ^ 2 '■ It ' ls obvious from Fig. 2a that ~ L x /\\\ magnetic field lines will 
undergo reconnection simultaneously (compared to a one by one line reconnection process 
for the Sweet-Parker scheme). Therefore the overall reconnection rate may be as large as 
V re c,giobai ~ Va (L x / Ay ) ( V^A \\ / rj) ^Z 2 , which means that the reconnection efficiency critically 
depends on the value of Ay. More realistically we will find that there are other global 
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constraints which end up determining the actual global reconnection speed. 

The relevant values of Ay and (y 2 ) 1 ^ 2 depend critically on the magnetic field statistics. 
Therefore in the next section we will briefly explore the expected properties of magnetic 
turbulence. 

3. Magnetic Field Structure in a Turbulent Medium 

The study of MHD turbulence is hampered by the absence of easily observable 
examples. We will base our discussion here on a model first proposed by Goldreich and 
Sridhar (1995, hereinafter GS95), which is at least consistent with available observations of 
the interstellar medium and the solar wind. A discussion of this model, and the justification 
for using it as a first approximation is given in Appendix A. In any case, the qualitative 
nature of our results, that a weak stochastic component to the field structure can have 
a dramatic effect on reconnection rates, is not sensitive to the details of the model we 
adopt. We discuss reconnection for various power spectra including the Kraichnan spectrum 
(Iroshnikov 1963, Kraichnan 1965) in Appendix D. 

For strong turbulence we have Aj_ ~ A||0(A||), although we can also define a weakly 
turbulent regime where the correlation length is much greater than the distance by 
which individual field lines deviate from a straight line. In either case, we can define a 
corresponding perpendicular wavenumber, kj_ ~ AT^Ay). 

In order to describe this wandering, we need to know something about quasi-static 
distortions of the field lines, or 'zero-frequency' modes. Normally, this means only those 
distortions of the magnetic field whose evolution is driven primarily by nonlinear effects. 
Here, however, we mean simply all modes that evolve on time scales long enough that 
they can be considered as a background to the evolution of flows driven by reconnection. 
Motions that can be described as Alfven waves may be relevant in this sense, as long as 
they involve very long wavelengths parallel to the magnetic field. However, most of our 
discussion will be concerned with modes with nonlinear time scales comparable to their 
wave periods, so that they are Alfven waves only in a very general sense. For simplicity, 
we will consider incompressible turbulence and assume that the magnetic field filling factor 
is not too small. The former should be a reasonable approximation for small scale, and 
therefore weak, fluid motions. The latter is probably a bad approximation inside the Sun, 
although reasonable whenever the magnetic pressure is comparable to, or greater than, the 
gas pressure (or when intermittency has not had time to develop). 

We can describe a model of MHD turbulence in terms of the relationship between 
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k± and hi at every scale in the turbulent cascade, as well as the scaling of as function 
of wavenumber. By restricting ourselves to a model in which field stochasticity is driven 
locally by a turbulent cascade we are neglecting the possibility that the field lines are 
anchored in some denser medium, like field lines extruded from the photosphere of the Sun, 
but even in this case the field lines are mixed at their footpoints by turbulent motions in 
the photosphere (see, for example, Parker 1988, and references contained therein). More 
generally there will be cases in which the level of stochasticity will be determined by 
nonlocal processes. Here we will restrict ourselves to cases where the field stochasticity can 
be understood in terms of local turbulence. 



If energy is injected on some scale /, with vi < Va, then GS95 predict that 

'vi_ 



*„ « r\kj)><* ( £ V /3 , (l) 



while the rms fluid velocity is given by 



«* « «i(*±0" 1/8 (^"J • (3) 

These equations are derived in Appendix A. Equation ([]]) is particularly interesting. It 
gives the geometry of small scale magnetic field structure in a turbulent cascade. This turns 
out to be critically important in magnetic reconnection. 

A bundle of field lines confined within a region of width y at some particular point 
will spread out perpendicular to the mean magnetic field direction as one moves in either 
direction parallel to the local magnetic field lines. The rate of field line diffusion is given 
approximately by 

d{y 2 } (y 2 ) 



dx An 



(4) 



where Ajj" 1 m k\\, k\\ is the parallel wavevector chosen so that the corresponding vertical 
wavelength, k±(k±), is ~ (l/ 2 ) 1//2 , and x is the distance along an axis parallel to the mean 
magnetic field. Therefore, using equation (|I|) one gets 

W ,Ml"f^ (5 ) 



dx \ I 2 ) \V A J 

where we have substituted (y 2 )^ 1 ^ 2 for k±. This expression for the diffusion coefficient will 
only apply when y is small enough for us to use the strong turbulence scaling relations, 
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or in other words when (y 2 ) < 1 2 (vi/Va) . Larger bundles will diffuse at a maximum 
rate of 1(vi/Va) 4 - For (y 2 ) small equation (|5|) implies that a given field line will wander 
perpendicular to the mean field line direction by an average amount 



I 1 ' 2 \V A . 

in a distance x. The fact that the rms perpendicular displacement grows faster than x is 
significant. It implies that if we consider a reconnection zone, a given magnetic flux element 
that wanders out of the zone has only a small probability of wandering back into it. 



4. Turbulent Reconnection 

4.1. Constraints on reconnection rate 

Outflow of matter from the reconnection layer constrains the achievable reconnection 
rates. In the presence of turbulence the thickness of the outflow layer increases with L x 
according to equation @: 

L r \ l l 2 ( vi \ 2 



when I > L x and 

when L x > I. Therefore the upper limit on V rec imposed by large scale field line diffusion is 



V rec < Va min 




Vl\' 

V A ) 



(9) 



This limit on the reconnection speed is fast, both in the sense that it does not depend 
on the resistivity, and in the sense that it represents a large fraction of the Alfven speed. 
Whether or not this limit is attainable depends to a large extent on the motion of the shared 
component of the magnetic flux. Nevertheless, we shall see below that equation (|^) is not 
merely an upper limit on the global reconnection speed, but under most circumstances, 
constitutes a reasonable estimate for its actual value. 

There are two important issues that we have neglected in arriving at this estimate. 
First, we have assumed an outflow velocity of Va spread over the entire width of the outflow 
region. However, in general the field lines that pass through the reconnection zone at any 
one moment will not fill the entire width of the outflow region. Instead, they will fill some 
fraction of it, with a proportionate drop in the effective outflow velocity. This reduction 



-10- 



will be almost exactly offset by the fact that every flux element stays in the reconnection 
zone for roughly L x /Va and that defines the meaning of "simultaneous" in this case. When 
reconnection starts, the outflow fills only a small fraction of the reconnection layer and 
this corresponds to the minimal estimate of the reconnection rate below. However, as 
reconnection proceeds, more independent flux elements become involved and this enhances 
the rate. 

In addition, equation (|9|) is based on the assumption that the field line topology far 
from the reconnection layer is not affected by the process of reconnection. In reality we 
expect that, once reconnection begins, the changing field line stresses will be communicated 
throughout the volume at a speed Va- Also, the reconnection layer will act as a source of 
small scale nonlinear Alfven waves. The outcome of these effects is somewhat uncertain, 
but it will probably have the effect of accelerating reconnection. 

We can obtain a lower limit on the reconnection speed by narrowing our focus and 
remembering that plasma can escape from the reconnection layer proper only by traveling a 
distance Ay (k±) along the local field lines, where fcj 1 ~ A is the thickness of the reconnection 
layer. A lower limit on the reconnection speed follows from the mass conservation condition: 

Vrec.local 

Ay « V A A, (10) 

and the requirement that at a fundamental level reconnection is driven by Ohmic diffusion: 

^rec local Ktri/A. (11) 

We have denoted this estimate for the reconnection speed with a special subscript not only 
because its value is determined by local physics, but also because we shall see that it is 
likely that the global speed of reconnection through a field line bundle is different, and 
much faster, than this local value. Combining equations ([!]), ([10]) and flTT|) and recalling 
that kj 1 ~ A and k^ 1 ~ Ay we get reconnection at a speed 



which implies that 

In arriving at this estimate we have treated field perturbations with frequencies less than 



V re c,local ~ V t . (13) 



~ V^/ An as essentially static. However, since we see from equations (10) and (11) that the 
reconnection rate r]/A 2 is ~ V^/Aii, this is a reasonable approximation. Finally, we note 
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that this process is controlled by field fluctuations with a perpendicular wavenumber 

, Vrec local ,_] '^l 

k± « ; ~ I — — . 14 

V Va\ v J 



We see from equation ([All ) that this is the largest wavenumber we can expect in the 



turbulent cascade, or at least the largest we can expect before the reconnection process 
generates its own local turbulence. 

If, however, the medium is partially ionized the maximum wavenumber given in 
equation ( |A11| ) should be reduced by a factor ~ (v/Vamb) 3 ^- This will reduce the 



reconnection speed by a factor of ~ (v/Vamb) 1 ^, but this reduction will be somewhat offset 
by the fact that the ions will be compressed until their pressure alone is comparable to the 
magnetic field pressure, and it is this compressed layer of ions which is ejected from the 
reconnection region along the field lines. The resulting velocity should be compared with 
the ambipolar reconnection velocity from (Vishniac & Lazarian 1998): 

/ 2r? \ 1/2 1 

VamUyolar = ~ Wl + 2/3x) 1/2 (15) 

V 7~recomb ' 1\ 

whether X is the ionization ratio, j3 is the ratio of the gas pressure to the magnetic 
pressure. The recombination time of the ions T recom b as well as Ohmic diffusivity r\ are 
taken far gas from the reconnection layer. The larger of the two velocities will determine 
the reconnection speed. Results in Table 1 in Vishniac & Lazarian (1998) show that the 
reconnection speed can be increased via ambipolar diffusion by a factor of 10 3 , which is not 
much, if we account for the fact that the original Sweet-Parker rate is very slow. A detailed 
account of reconnection in the presence of MHD turbulence and ambipolar diffusion is given 
in Lazarian & Vishniac (1998). 

So far we have ignored the fact that the thin current sheets envisioned in Sweet-Parker 
reconnection are unstable. We briefly address this issue in Appendix C and show that 
allowing for the presence of tearing modes raises the minimum value of V rec ^ oca i to 

3/4 3/16 

~ V A (j^J (y^l) • Even in the context of large scale turbulence in the interstellar 
medium, this is only an increase of about one order of magnitude, and still leaves us with 

<^v A . 



4.2. Estimating the global value of V r 



The physics of the reconnection layer will be significantly more complicated than the 
simple picture above indicates. First of all, a number of known processes can enhance 
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reconnection. In the appendices we discuss Bohm diffusion, anomalous resistivity and tearing 
modes. All of these have been discussed previously as ways to accelerate reconnection, but 
none have emerged as universal mechanisms for producing fast reconnection. In a previous 
paper (Vishniac and Lazarian 1998) we have shown that while ambipolar diffusion can 
produce dramatically faster reconnection rates, the overall pace of reconnection is still orders 
of magnitude below Va- Anomalous resistivity (see Appendix B) is of marginal importance 
for interstellar reconnection, although it can be important for laboratory plasma. Finally, 
Bohm diffusion (Appendix B) can provide a strong enhancement of reconnection speeds, 
but unfortunately it is unclear whether or not the empirical concept of Bohm diffusion is 
applicable to astrophysical plasmas. 

All of this is probably less important than the global ordering of reconnection events. 
The minimal estimate of V rec given in the previous subsection is based on the assumption 
that reconnection proceeds sequentially, that is, the reconnection speed is simply the speed 
with which reconnection propagates through a single flux element. This is not obviously 
correct, since the reconnection zone contains many independent reconnection events at any 
one time. We need to define a global reconnection speed, V reC)g i b a i, which describes the 
rate at which flux is reconnected throughout the reconnection zone. In order to arrive at 
a reasonable estimate of this speed, we have to determine which aspect of reconnection 
sets a limit on its efficiency. Based on our previous discussion, we can see that there are 
four possibilities: the mass flow from the reconnection zone itself, the speed with which 
reconnected flux elements move across the reconnection zone and off the edge, the ejection 
of the flux associated with the shared magnetic field component, and the mass flow from 
the contact volume (roughly everything within a distance L x of the reconnection zone). In 
the case of Sweet-Parker reconnection the first process provides the critical constraint (and 
the third and fourth are not separate constraints). 

The first limit is determined by V reCt i oca i, but is not equal to it. Since the current sheet 
contains as many as k±k\\L x independent reconnection surfaces, each one reconnecting flux 
at a rate \± V reC) i oca i, the global reconnection speed could be enhanced by a factoif|of kuL x . 
In other words, using equations (H), fll3D, and ([14]) we get 



8 The ambipolar reconnection rate given by Eq. ( p~5| ) is a local reconnection speed. Therefore, in the 
presence of turbulence the global ambipolar reconnection speed is larger by a similar factor (Lazarian & 
Vishniac 1999). 





Since this will almost always be less restrictive than equation (|9|), we conclude that the flow 
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of matter from the reconnection zone itself does not limit the speed of reconnection. 

The second limit is set by the speed with which reconnected magnetic field elements 
can pass through one another on their way to the edge of the reconnection zone. After 
initial reconnection magnetic field elements will move until they encounter magnetic flux 
moving in the opposite direction. If this process limits the speed of reconnection, then 
we can assume that each flux element has a typical size comparable to the width of the 
reconnection zone, and moves a similar distance before being entangled by another flux 
element. (If the reconnected flux elements were distributed more sparsely, then from the 
preceding paragraph we can see that the gaps would fill in with freshly reconnected flux.) 
We can restrict ourselves to the case where the flux elements cross the reconnection layer 
vertically, rather than at a shallow angle, since in the latter case reconnection between 
opposing flux surfaces would, if anything, be enhanced and the volume filled by a single 
reconnected element would evolve into one filled by many flux elements with a nearly 
vertical orientation. Outside the reconnection layer the magnetic field moves parallel to the 
reconnection layer only a small amount, and reverses direction after each interchange. The 
layers of magnetic field around the reconnection zone are removed starting from its edges 
and working inwards. 

The steady pace of reconnection ensures that the reconnection zone is in a constant 
state of turbulent motion. In the absence of collisions between flux elements we could 
assume an average transverse velocity of Va, but in their presence we need to estimate 
the reconnection speed for these flux elements. In an important sense, this is just a 
repetition of the problem we face in finding V reCtg i b a i- The only difference is that turbulent 
motions within the reconnection zone guarantee a local source of turbulence. The level of 
turbulence in a reconnection layer of width A is given by equating the energy injection rate, 
Vjl(V e j ect / A), with the turbulent dissipation rate, ~ tf^/V^A. We conclude that 



We note that the limiting value of k± imposed by dissipation outside the reconnection zone 
will be replaced locally by some much larger value due to this separate source of turbulence. 
Consequently, if the motion of flux elements in the reconnection zone limit the global pace 



^local 




(17) 



of reconnection, and if we adopt the ansatz^j that 




(18) 



9 We do not consider the possibility of a logarithmic dependence on the Lundquist number, as in this case 
reconnection is already fast. 
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where m and n are to be determined, then the average ejection speed of magnetic flux from 
the sides of the reconnect ion zone is just 

^MM^vrV- (19) 

On the other hand, if reconnected magnetic flux elements fill a large fraction of the 
reconnection zone, as they must if they determine the overall pace of reconnection, then 

Vrec,global ~ Reject- (^0) 

We see from equations and (p0|) that 



V recSohal « V A (^rj " (m/4) • (21) 

We can see from this that m < 4, but there is no obvious constraint on n. However, since 
the intersection of separate flux elements constitutes a repetition, on a smaller scale, of the 
global process of reconnection, we expect to find a still smaller scale A' which characterizes 
the size of the flux elements generated during the reconnection of flux elements of size A. 
The same arguments then allow us to conclude that 



ri \ l-(m/4) 

V rec , global « V A { -^j . (22) 



From equations fl2l"|) and fl22|) we see that n = 0, and (now referring to equations (|T 
and (^)) that m = as well. In other words, adopting the hypothesis that reconnection 
is limited by rate of interchange of the tangle of reconnected flux elements that cross the 
reconnection layer we conclude only that V reCt giobai ^ Va- Since this limit is less stringent 
than equation (|9|) we see that this process is not, in fact, the bottleneck for turbulent 
reconnection and may therefore be ignored. 



,max ' 



Our argument assumes that the scale A is greater than the spectrum cutoff, fcj 1 ^ 
for the external turbulence and that A' is similarly greater than the local spectrum cutoff 
k' ± max . We can check this by adopting the very conservative constraint that the flux element 
widths are greater than the scale imposed by Ohmic diffusion with the time necessary for a 
reconnection front to propagate over the reconnection layer, that is 

A > ( " L - ) U \ V2 L. ( VA ) , (23) 

\ ^rec, global 

I \L x VaJ \v rec 

,global / 



and 



A'>MM" 2 , (24) 



V., 



rec, global / 
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or 

A' / y \ V2 

A V AV rec ,global / 

From equation ([All ) we have 

,„ 4> f*)'"f ^'"(M"^^"* (26) 



T] J \ ^rec, global J \ ^rec, global J \ I 

Since all these factors, except possibly the last, are larger than one, and the first is usually 
much larger than one, we can conclude that k^^ max A > 1. Similarly, the local turbulence 
gives us a modified limit on k± which is 

U A _i f Vre^globa^ 1 ^ ( Va^ \ ^ 

Combining this with equation (|2"5]) we have 



*- a->( i )' / 7M) 3/ Vw^V /4 > fM) 1/ 7^^) 1/4 (28) 

,global J \ V 

J \ V A J \ V J \Vr ■ec, global J 

which is strictly greater than one. We conclude that the our assumption of a self-similar 
hierarchy is self-consistent. 

We note that the fact that the horizontal transport of flux elements does not provide 
the limiting constraint on V r ec,giobai implies that the fraction of the reconnection zone filled 
with already reconnected flux elements, f r , is not necessarily large. However, there is firm 
lower limit on f r . Since V rec>g i obal = f r V eject then since V eject < V A , 

r ^ Vrec,global / 9q n 
Jr> y A ■ {^) 

If reconnection is fast then f r must be of order unity. 

The third limit involves the width of the zone through which the shared component 
of the magnetic field is ejected, w e j ect . Given the stochastic geometry of the initial field 
configuration, and its evolution during reconnection, this is not necessarily the same as the 
width of the reconnection zone. In any case we have 

T/ T _ T/ _ Vrec,global fon\ 

Vrec^lobal-I^x ~ V eject™ eject ~ W , {o(J) 

Jr 

or 

W e ject — frL x . (31) 
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In a certain sense, this is just a restatement of the hypothesis that the ejection of the shared 
flux from the contact volume controls the speed of reconnection. However, there is another 
connection between w e j ect and f r where causality runs in the other direction. If we consider 
the initial surface layer of flux elements facing the reconnection zone, they have connections 
which run deep into the underlying magnetized region (cf. equations (0) and (fD). However, 
these connections have a small volume filling factor. Once these initial elements reconnect, 
further reconnection involves flux elements that are pressed down between the reconnected 
elements. This gives them an extra component normal to the plane of the reconnection 
zone, which we can express as a displacement of the flux element from its initial position. 
As f r — > 1 (or at least a fraction of order unity) this displacement is typically as large 
as the diffusion limit given by equations (|7|) and (0) and the volume filling factor of the 
reconnected flux elements rises to one. Since typical flux elements will acquire a vertical 
stretch equal to this displacement, we see that this will also be the width of the escape zone 
for the shared flux. In other words, if reconnected flux elements fill up the reconnection 
zone, then w e j ect will be of order f r and shared flux can escape without any difficulty 

Is this a likely scenario? If the reconnection zone is sparsely filled, then nearby flux 
elements will be forced into it by the local pressure excess. Since matter can escape without 
difficulty in this model, this pressure excess is equal to the magnetic pressure. This will 
not be sufficient to force magnetic field lines to bend at sharp angles, unless the whole 
excess is supported on a few field lines. However, it will be sufficient to force nearly parallel 
field lines to move, so we can envision this motion as consisting almost entirely of motion 
towards the reconnection zone, without any significant bending away from the initial mean 
field direction parallel to the reconnection zone. Furthermore, as long as w e j ec t is no larger 
than the value of /cj 1 corresponding to parallel wavelengths ~ L x , then the turbulent 
motions within the magnetized fluid can be assumed to keep moving flux elements up to 
the reconnection zone. This will ensure that the reconnection zone is continuously supplied 
with magnetic elements to reconnect, and consequently the f r will grow until the loss of 
magnetic flux from the reconnection zone balances the supply. This implies that w e j ec t will 
grow until some other process limits the global speed of reconnection. 

The fourth, and last, limit is equation (|9|). We see from the three preceding arguments 
that this is not just an upper limit, but a reasonable estimate for the global speed of 
reconnection in a turbulent medium. Reconnection will generally occur at a substantial 
fraction of the Alfven speed with the level of turbulence in the medium (or alternatively, 
the level of stochasticity in the magnetic field) controlling the exact speed. 

We can turn this line of reasoning around, and use the global speed of reconnection to 
estimate the properties of the reconnection zone itself. First, we note that the field lines 
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will bend down to the reconnection zone with a typical displacement equal to the diffusion 
distance, or for L x < I, 



Vl 



2 £3/2 



W -** = W w" (32) 

This represents a systematic bending which is just big enough to compete with random 
diffusion near the central parts of the contact volume. Near the edges of the reconnection 
zone it will be much more conspicuous. Also by our previous reasoning, we can expect that 

A~(£Vfer<i- (33) 



v A J V / 

The field lines will bend drastically, with a large concentration of current and Ohmic 
dissipation, only in the reconnection zone itself. If we assume that the reconnection zone 
and current sheet are identical then 

Vrec,local = (34) 

This is probably not quite right, but the effect of tearing modes appears to be only a 
modest enhancement of reconnection rates, judging from the results of Appendix C, and 
the turbulent motions of reconnected field lines seem likely to disrupt modes with k\\ <C k±. 
We should nevertheless bear in mind that equation ([34]) may lead to an underestimate of 
the reconnection zone width, and an overestimate of Ohmic heating within it. 

Since global reconnection speeds are given by equation (Q), and since the accumulation 
of reconnected flux elements in the reconnection zone does not limit the pace of reconnection, 
we can use 

V re c,global ~ k\\LxVrec,local ~ k±k\\L x T), (35) 

where, as before, we have written the width of the reconnection zone in terms of the 
equivalent wavenumber in the turbulent cascade outside the reconnection zone. Using 
equation ([!]) we can rewrite this expression and equate it to the global reconnection speed 
given in equation We find 



^'{tJ UfJ If J • < 36 > 

where we have assumed that L x < I, for simplicity, and because this seems like the most 
likely case. The Ohmic dissipation rate, f]/A 2 , is actually less than kvVj^ because the 
pace of reconnection is set by the escape of matter from the contact volume, rather than 
the reconnection zone. Consequently the ejection of matter from the reconnection zone 
is inhibited and the corresponding ejection velocity is less than Va- In other words, the 
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build-up of excess matter throughout the contact volume softens the pressure gradient 
so that escaping fluid reaches an ejection velocity ~ Va only after it has already left the 
reconnection zone. 

One of the more striking aspects of our result of the global reconnection speed is that it 
is relatively insensitive to the actual physics of reconnection. Equation (§) only depends on 
the nature of the turbulent cascade. Although we have reached this conclusion by invoking 
a particular model for the strong turbulent cascade, we can see that any sensible model will 
give qualitatively similar results. 

The conclusion that reconnection is fast, even when the local reconnection speed is 
slow, represents a triumph of global geometry over the slow pace of Ohmic diffusion. In the 
end, reconnection can be fast because if we consider any particular flux element inside the 
contact volume, assumed to be of order L\, the fraction of the flux element that actually 
undergoes microscopic reconnection vanishes as the resistivity goes to zero. 

4.3. Reconnection for a modified spectrum of turbulence 

Our main conclusion, that reconnection is dramatically enhanced in the presence of 
a stochastic field, does not depend on the particular model of turbulence we adopt. Our 
work here assumes a locally generated turbulent spectrum which follows the scaling laws 
suggested by GS95. However, we depend only on the scaling law for power as a function 
of k± and the eddy anisotropy, that is k»(k±). The former is supported by observations of 
the solar wind. The latter depends on the assumption of locality of nonlinear interactions 
in wavenumber space and may be wrong (cf. Matthaeus et al. 1998). In Appendix D we 
discuss magnetic reconnection for an arbitrary power-law relation between k\\ and k± and 
conclude that if fcn ~ k p ± with p > 1/2 then the reconnection rate is independent of the 
Lundquist number. 

Consider now the processes that can modify the spectrum of turbulence. The turbulent 
spectrum may be changed due to: (a) a weak large scale magnetic field, (b) the injection 
of energy on small scales as a result of reconnection, (c) large scale discontinuities resulting 
from bringing magnetic fields with different directions into direct contact (by definition this 
is a necessary part of reconnection), (d) the backreaction of enhanced reconnection on the 
inverse cascade, and (e) the injection of field stochasticity from some other region. 

We discussed (a) and (b) above. The first simply makes the field more chaotic and 
enhances reconnection rates, since the effective size of I is reduced. The latter has a much 
weaker effect, since it has only a slight effect on the large scale field structure, but is likely to 
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marginally increase reconnection speeds. Apropos (c), we expect that for strong magnetic 
fields brought into contact the efficiency of the coupling between eddies at different sides 
of the reconnection layer will be low. Consider, for example, two flux tubes at 7r/2 angle. 
The eddies with k\\ }Up in the upper flux tube will be exciting perturbations of the order fcjj~* p 
perpendicular to the magnetic field lines in the lower flux tube. However the time scales 
for perturbations with kx,i<m> — k\\ up will be very different from the time scale of the eddies 
providing excitation. Due to this mismatch in the time scales the crosstalk between eddies 
on opposite sides of the reconnection zone will be suppressed. 

The high speed of reconnection given by equation (|D naturally leads to a question 
of self-consistency (our point (d) above). Is it reasonable to take the turbulent cascade 
suggested in GS95 when field lines in adjacent eddies are capable of reconnecting? It turns 
out that in this context, our estimate for V re c,giobai is just fast enough to be interesting. We 
note that when considering the intersection of nearly parallel field lines in adjacent eddies 
the acceleration of plasma from the reconnection layer due to the pressure gradient is not 
k\\V%, but rather OfAD^j, since only the energy of the component of the magnetic field 
which is not shared is available to drive the outflow. On the other hand, the characteristic 
length contraction of a given field line due to reconnection between adjacent eddies is only 
k\\jk\. This gives an effective ejection rate of fen Va. Since the width of the diffusion layer 
over a length A;^" 1 is just , we can replace equation (Q) with 

k\\ 

Vrec,global ~ V AT~ ■ (37) 
K± 

The associated reconnection rate is just 

T relonnect ~ V Ak\\, (38) 

which in GS95 is just the nonlinear cascade rate on the scale k7 . However, this result is 
general and does not involve assuming that GS95 is correct. We will see in §5.1 that most 
of the energy liberated in reconnection goes into motions on length scales comparable to 
the dimensions of the reconnecting eddies, so this energy release will not short circuit the 
energy cascade described in GS95. On the other hand, we can invert this argument to see 
that reconnection can play an important role in preventing the buildup of unresolved knots 
in the magnetic field. Such structures could play a major role in inhibiting the cascade of 
energy to smaller scales, flattening the energy spectrum relative to the predictions of GS95. 
Our conclusion is that such structures will disappear as fast as they appear, supporting the 
notion that they play a limited role in the dynamics of MHD turbulence. 

Finally, we note that if the magnetic field structure is driven by turbulence in another 
location, as when the footpoints of magnetic arcades are stirred by turbulent motions, then 
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we can evaluate its effects in terms of the amplitude of field stochasticity and the scaling 
of structure anisotropy with scale. In the absence of any particular model for this process, 
we note that the robust nature of our conclusions, as shown in Appendix D, implies that 
under these circumstances reconnection will be sensitive to the amplitude of the induced 
field stochasticity, but not the details of the turbulent mixing process. 



5. Discussion 

5.1. Energy dissipation 

The usual assumption for energy dissipation in reconnection is that some large fraction 
of the energy given up by the magnetic field, in this case ~ pV\L\^ goes into heating the 
electrons. This is not the case here. Only a fraction, ~ l/(k\\L x ) of any flux element is 
annihilated by Ohmic heating within the reconnection zone. Over the entire course of the 
reconnection event the efficiency for electron heating is no greater than 

77/ A V rec 

,local J- / on \ 

e e ^ 77 = TP = -j—j- ■ l dy J 

'rec, global 'rec, global ™\\*-'x 

or, from equations (H), (|l|), and (|36|) , 

^rvM^'f . (40) 



V J W A J 

The electron heating within the current sheet will not be uniform, due to the presence of 
turbulence, the intermittent presence of reconnected flux, and any collective effects we have 
neglected here. To the extent that these are important they will also lower the electron 
heating efficiency by broadening the reconnection layer. 

Naturally, the low value of e e is of little interest when ion and electron temperatures 
are tightly coupled. When this is not the case our model for reconnection may have 
testable consequences. As an example, we may consider advective accretion flows (ADAFs), 
following the general description given in Narayan and Yi (1995) in which advective flows 
can be geometrically thick, and optically thin with a small fraction of the dissipation going 
into electron heating. If, as expected, the magnetic pressure is comparable to the gas 
pressure in these systems, then a large fraction of the orbital energy dissipation occurs 
through reconnection events. If a large fraction of this energy goes into electron heating (cf. 
Bisnovatyi-Kogan fc Lovelace 199"T|) then the observational arguments in favor of ADAFs 



are largely invalidated. Our results suggest that reconnection, by itself, will not result in 



-21 - 



channeling more than a small fraction of the energy into electron heating 

We can also estimate the efficiency of the conversion of magnetic field energy to high 
frequency Alfven waves. One channel for the creation of such waves will be as waves 
radiated from the turbulence in the reconnection zone. The rate of energy emission is hard 
to estimate since the nature of the turbulence in the reconnection zone is so uncertain. 
Plausibly the energy lost from the this zone in the form of radiated Alfven waves will be 
at least comparable to the Ohmic dissipation, but it is difficult to make a more precise 
estimate. On the other hand, the region immediately outside the reconnection zone will 
also be affected, via the constant interchange of reconnected flux and the subsequent 
relaxation of field lines just outside the reconnection zone. The total energy involved at 
any one time will be ~ pV\L\!S. and the rate of excitation of Alfven modes should be of 
order k\\V e j ect = k\\V reCt giobai- This process will operate for a time ~ L x /V reCt gi b a i, so that the 
efficiency of wave energy generation will be 

(_n_ V' 5 (Va\ 2 (^} 1,m 



- ^ Ak « ~ iti) lifJ W ■ (41) 

The associated wave frequencies will be distributed between k\\V reC! gi b a i and Va/A depending 
on the nature of the turbulent zone, that is 

V A (V A lV /s ( v, \*/» / 1 \ 1/5 V A (VaIV" I v, \W ( I \ 3/1 " , M 

These limits differ by a factor of (VaI/t]) 1 ^ 5 , which is large enough to make precise 
predictions difficult, but small enough that we have a qualitative sense of the results. 

These waves are unlikely to affect of the topology of the magnetic field, since their 
typical scale is, at most, only Ay [_]• However, they may play a role in the anisotropic heating 
of ions in stellar coronae. In this connection we note recent observations of anisotropic 
velocity distributions for HII and OVI ions in the solar corona ( [Kohl et al. 1998| ) which they 
suggest may be due to the presence of high frequency (tens of kilohertz) Alfven waves (see 



Cranmer, Field, fc Kohl 1998 for a detailed analysis). 



10 Othcr authors (e.g. Gruzinov 1998) have suggested that MHD turbulence in a collisionless, highly 
magnetized environment may not terminate on small scales in ion heating, so that ADAFs will necessarily 
deposit a large fraction of their heat into electrons. However, Gruzinov's calculation rests on treating the 
turbulence as a collection of independent linear waves, which is an unconvincing model for disturbances with 
a coherence time comparable to a single wave period. 

11 However, the large scale relaxation of magnetic field lines is likely to feed energy back into the turbulent 
cascade, providing a nonlinear feedback mechanism which will accelerate reconnection when the initial field 
stochasticity is small. 



-22 - 



5.2. Dynamos 

To enable sustainable dynamo action and, for example, generate a galactic magnetic 
field, it is necessary to reconnect and rearrange magnetic flux on a scale similar to a galactic 
disc thickness within roughly a galactic turnover time (~ 10 8 years). This implies that 
reconnection must occur at a substantial fraction of the Alfven velocity. The preceding 
arguments indicate that such reconnection velocities should be attainable if we allow for a 
realistic magnetic field structure, one that includes both random and regular fields. 

One of the arguments against traditional mean-field dynamo theory is that the rapid 
generation of small scale magnetic fields suppresses further dynamo action (e.g. Kulsrud 
& Anderson 1992). Our results thus far show that a random magnetic field enhances 
reconnection by enabling more efficient diffusion of matter from the reconnection layer. 
This suggests that the existence of small scale magnetic turbulence is a prerequisite for 
a successful large scale dynamo. In other words, we are arguing for the existence of 
a kind of negative feed-back. If the magnetic field is too smooth, reconnection speeds 
decrease and the field becomes more tangled. If the field is extremely chaotic, reconnection 
speeds increase, making the field smoother. We note that it is common knowledge that 
magnetic reconnection can sometimes be quick and sometimes be slow. For instance, the 
existence bundles of flux tubes of opposite polarity in the solar convection zone indicates 
that reconnection can be very slow. At the same time, solar flaring suggests very rapid 
reconnection rates. 

There is a general belief that magnetic dynamos operate in stars, galaxies (see Parker 
1979) and accretion disks (see Balbus & Hawley 1998). In stars, and in many accretion 
disks, the plasma has a high ft, that is the average plasma pressure is higher than the 
average magnetic pressure. In such situations the high diffusivity of magnetic field can 
be explained by concentrating flux in tubes (Vishniac 1995a,b). This trick does not work 
in the disks of galaxies, where the magnetic field is mostly diffuse (compare Subramanian 
1998) and ambipolar diffusion impedes the formation of flux tubes (Lazarian & Vishniac 
1996). This is the situation where our current treatment of magnetic reconnection is most 
relevant. However, our results suggest that magnetic reconnection proceeds regardless 
and that the concentration of magnetic flux in flux tubes via turbulent pumping is not a 
necessary requirement for successful dynamos stars and accretion discs. 

We do not address here the controversial issue of the turbulent dynamo in clusters 
of galaxies. This was first suggested by (Jaffe 1980) and was elaborated in great detail 
by Ruzmaikin, Sokoloff & Shukurov (1989), who claimed an excellent match between 
observations and predictions based on the Kazantsev (1968) theory of the turbulent 
dynamo. However, Goldshmidt & Rephaeli (1993) found a large (~ 10 20 ) numerical error 
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in the value of Ohmic diffusivity used by Ruzmaikin, Sokoloff & Shukurov (1989), which 
formally invalidated their result. However, if it were possible to use the effective diffusivity 
determined by the reconnection rate instead of Ohmic diffusivity, then the theory of 
turbulent dynamo can be revived for clusters of galaxies. We postpone the discussion of 
this important question until our next paper. 

Our results show that in the presence of MHD turbulence magnetic reconnection 
is fast, and this in turn allows the possibility of 'fast' dynamos in astrophysics (see the 
discussion of the fast dynamo in Parker (1992)). As the fluid conductivity tends to zero the 
relevant parallel wavelength also goes to zero, allowing a change of magnetic field topology 
via reconnection at a finite rate. 



5.3. Turbulent diffusion, particle diffusion, and self-excited reconnection 

It may seem from our discussion that this paper revives the concept of turbulent 
diffusion. However, a closer comparison between the concept original concept of turbulent 
diffusivity and our suggestion shows that the two are very different. Within the turbulent 
diffusivity paradigm, magnetic fields of different polarity were believed to filament and 
intermix on very small scales while reconnection was believed to proceed slowly. On the 
contrary, we have shown that the global speed of reconnection is fast if a moderate degree 
of magnetic field line wandering is allowed. The latter, unlike the former, corresponds to a 
realistic picture of MHD turbulence and does not entail prohibitively high magnetic field 
energies at small scales. 

On the other hand, the diffusion of particles through a magnetized plasma is greatly 
enhanced when the field is mildly stochastic. There is an analogy between the reconnection 
problem and the diffusion of cosmic rays (Barghouty & Jokipii 1996). In both cases charged 
particles follow magnetic field lines and in both cases the wandering of the magnetic field 
lines leads to efficient diffusion. 

Finally, we note that observations of Solar flaring seem to show that such events 
start from some limited volume and spread as though a chain reaction from the initial 
reconnection region initiated a dramatic change in the magnetic field properties. Indeed, 
Solar flaring happens as if the resistivity of plasma were increasing dramatically as plasma 
turbulence grows (see Dere 1996 and references therein). In our picture this is a consequence 
of the increased stochasticity of the field lines rather than any change in the local resistivity. 
The change in magnetic field topology that follows localized reconnection provides the 
energy necessary to feed a turbulent cascade in neighboring regions. This kind of nonlinear 
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feedback is also seen in the geomagnetic tail, where it has prompted the suggestion that 
reconnection is mediated by some kind of nonlinear instability built around modes with 
very small k\\ (cf. Chang 1998 and references therein). The most detailed exploration 
of nonlinear feedback can be found in the work of Matthaeus and Lamkin (1986), who 
showed that instabilities in narrow current sheets can sustain broadband turbulence in two 
dimensional simulations. Although our model is quite different, and relies on the three 
dimensional wandering of field lines to sustain fast reconnection, we note that the concept 
of a self-excited disturbance does carry over and may describe the evolution of reconnection 
between volumes with initially smooth magnetic fields. 

6. Conclusions 

In this paper we have explored the consequences of a small stochastic component for 
the reconnection of large scale magnetic fields. We have used a particular model of MHD 
turbulence, due to Goldreich and Sridhar (1995), in order to draw firm conclusions, but 
the main features of our results should follow from any model which is consistent with our 
observations of MHD turbulence in astrophysical systems. We have implicitly assumed that 
viscosity is less important than resistivity in setting the limits for any turbulent cascade. 
Although this should usually be the case for astrophysical plasmas, it is not guaranteed, and 
the effects of viscosity will need to be considered in more detail (Lazarian & Vishniac 1999). 
Finally, we have assumed that we are dealing with a strong magnetic field, where motions 
that tend to mix field lines of different orientations are largely suppressed. The galactic 
magnetic field is usually taken to have grown via dynamo action from some extremely 
weak seed field (cf. Zel'dovich, Ruzmaikin & Sokoloff 1983, Lazarian 1992, Kulsrud, Cen, 
Ostriker, & Ryu 1997 and references contained therein). When the field is weak it can 
be moved as a passive scalar and its spectrum will mimic that of Kolmogorov turbulence. 
The difference between A^ and An vanishes, the field becomes tangled on small scales, and 
Vrec,iocai becomes of the order of Va- Of course, in this stage of evolution Va may be very 
small. However on the small scales Va will grow to equipartition with turbulent velocities 
on the turnover time of the small eddies. The enhancement of reconnection as Va increases 
accelerates the inverse cascade as small magnetic loops merge to form larger ones. 

We list the principal results of this paper below. 

First, the rate of magnetic reconnection is increased dramatically in the presence of 
a stochastic component to the magnetic field. This component arises naturally whenever 
turbulence is present. Even when the turbulent cascade is weak the resulting reconnection 
speed is independent of the Ohmic resistivity. 
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Second, the argument that the rapid rise of random magnetic field associated with 
dynamo action results in the suppression of dynamo (Kulsrud & Anderson 1992) is 
untenable since the increase of the random component of the magnetic field increases the 
reconnection rate. We conclude that dynamo is a self-regulating process. 

Third, the second parameter in determining the reconnection speed is not some aspect 
of the microphysics, but the level of field stochasticity (or the large scale kinetic energy 
feeding the turbulent cascade). As reconnection proceeds the local turbulent cascade will 
grow stronger and the initial level of stochasticity will matter less and less. On the contrary, 
microphysical processes widely believed to speed up reconnection, i.e. anomalous resistivity, 
fail in interstellar conditions. 

Fourth, since the speed of reconnection depends on the presence of small scale 
turbulence, non-resistive processes which terminate the turbulent cascade may slow 
reconnection down. As an example, we can cite neutral-ion collisions in a dense, partially 
ionized gas. On the other hand, the rate of collisions in a fully ionized plasma is irrelevant. 
Under typical astrophysical conditions v <C r\ for fully ionized plasmas. When collisions 
become very rare, the only effect is to eliminate the cascade of compressive modes, which 
are not important for reconnection. 

Fifth, reconnection in a turbulent medium deposits only a small fraction of the 
magnetic field energy directly into heating electrons. Whether or not this is important 
depends on the rate of energy exchange between ions and electrons, and the unknown 
degree to which the turbulent energy cascade can avoid depositing its energy into electrons. 

Sixth, regardless of the other arguments advanced here, there is a minimal reconnection 
speed, ~ VaR<l 3 ^ W , where 1Z L = (VaI/t)), which is much faster than the Sweet-Parker 
estimate, even though it is still unrealistically slow. 

Finally, the current sheet associated with reconnection is not substantially broadened 
from the Sweet-Parker result, even when reconnection is fast. 
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A. The Power Spectrum of Magnetohydrodynamic Turbulence 

Here we discuss a model for MHD turbulence first proposed in GS95. Although their 
work included an explicit calculation of coupling constants based on treating the turbulent 
motions as shear Alfven waves, the important parts of the model depend only on some 
rather general scaling arguments, which we will repeat here. These arguments depend on 
assuming that power is transferred locally in Fourier space, which is a convenient, but not 
rigorously justified, assumption. 

We begin by noting that the first simple model of incompressible MHD turbulence 
was proposed independently by Iroshnikov (1963) and Kraichnan (1965) based on the 
interactions of triads of waves (see also Diamond and Craddock 1990). In this picture the 
energy transfer rate is roughly 

« ^ (AD 

where k is the magnitude of a wavevector, vj. is the rms fluid velocity contributed from 
power on the scale fc -1 , uja = k\\VA is the Alfven wave frequency, k\\ is the wavevector 
component parallel to the magnetic field direction, and Va is the Alfven velocity. In the 
original picture the power was assumed to spread isotropically in wavevector space. As long 
as v k < Va this expression for the nonlinear time scale will be less than kvk, that is, the 
magnetic field reduces the cascade of energy to higher wavenumbers. 

In a recent series of papers Goldreich and Sridhar (Sridhar and Goldreich 1994; 
Goldreich and Sridhar 1995; Goldreich and Sridhar 1997) questioned the validity of this 
picture and argued that the diffusion of power towards larger values of kn is strongly 
suppressed. Their original claim that three wave interactions are completely suppressed 
has been strongly criticized (Montgomery and Matthaeus 1995, Ng & Bhattacharjee 1996). 
However, in their most recent paper (Goldreich and Sridhar 1997) they showed convincingly 
that the effect of residual three wave couplings is consistent with a picture in which the 
basic nonlinear time scale is set by equation ( |A1~|) but with an anisotropic spectrum in 
which virtually all of the transfer of power between modes moves energy towards larger k± 
while leaving k\\ unchanged. They proposed calling this regime 'intermediate turbulence', 
since while the nonlinear decay rate is identical to the usual expression for weak turbulence 
among dispersive waves, in this case the higher order mode couplings are all comparably 
important. If we invoke the constancy of the local energy flux through the cascade, vl/r n i, 
as a function of scale, then from equation (|A1|) we see that in this regime 

v k oc kl l/ \ (A2) 

where we have assumed that k± » k\\. 
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As the power cascades to larger values of k± the magnetic field becomes progressively 
less important in the mode dynamics. Eventually we have 

k\\V A < kv k . (A3) 

In this limit the motions are no longer wave-like, the magnetic field exerts only a weak 
influence on the dynamics, and the fluid motions resemble ordinary hydrodynamical 
turbulence with a nonlinear time scale ~ kv k - Since feii is no longer privileged, the cascade 
of power is in the direction of increasing isotropy. Since k± 3> k\\ this implies an increase in 

fc||. 

If energy is injected on some scale I, with vi < Va, then we expect the cascade to transfer 
energy to larger k± until the condition expressed in equation ( |A3| ) is satisfied. At this point 
the turbulence is no longer wavelike (since r~l ~ ^Va). However, since turbulence tends 
towards isotropy when the magnetic field is completely negligible, we expect kn and k± to 
increase in tandem so that equation ([A3] ) is just marginally satisfied. This is the regime of 
strong turbulence described in GS95. At all smaller parallel wavelengths fluid motions bend 
magnetic field lines easily. Consequently we expect most of the power in energy spectrum 
to be centered around wavenumbers such that 

fell ~ k ± £- (A4) 
Va 

We approximate the energy transfer rate for the turbulent cascade, £, for v i < Va as 

i * (A5) 

if Vi < Va on the scale I. The usual hydrodynamic choice £ ~ -j- is valid otherwise, although 
not relevant for our present discussion. When the magnetic field is weak and the largest 
scales in the turbulent cascade are essentially hydrodynamic then we can identify I with the 
scale of equipartition so that vf ~ V\. 

Conservation of energy in the turbulent cascade implies that £ = v\jr n i = const and is 
only determined by large scale motions. If I is the large scale eddy size, and vi < Va then 
we use equation flA"5| ) to find £. From equation ( |A1| ) we find that 



IV A T n l k\\V A 



Combining equations ( |A^) and ( |Aq ) we obtain 

k^l-\kj)^[^f'\ (A7) 
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Substituting this in equation ( |A1|) gives the energy transfer rate 

1/3 



« ?(M 2/S (£) , (A8) 



while the rms fluid velocity on the scale k follows from equation flA4|) 

^~^(M" 1/3 (^-) 1/3 - (A9) 

Had we assumed an isotropic distribution of power in k space (that is, had we failed 
to distinguish between k± and k») we would have recovered the Kraichnan spectrum, with 
Vk oc £; _1//4 . Instead, we find v/~ oc k» Curiously enough, this anisotropic cascade looks 
qualitatively like the Kolmogorov spectrum when expressed in terms of k±. This becomes 
less surprising if one recalls that the magnetic field does not influence motions that do not 
bend it. We also note that for the characteristic k\\ given in equation ( |A7| ) the interaction 
between different modes can be viewed, with equal validity, as a wave-wave interaction, an 
eddy-eddy interaction, or a wave-eddy interaction. That the first two give roughly equal 
nonlinear time scales is immediately apparent from equations (|AT|), and ( |A4j) . That the 
third point of view is also equivalent follows from the work of Similon and Sudan (1989) 
who showed that the dissipation of Alfven waves in a stochastic field can be estimated by 
replacing the microscopic Ohmic resistivity by the diffusion coefficient for waves traveling 
along the field lines. In terms of the turbulent cascade envisioned here, this is equivalent to 
making the substitution 

V -> Vturb « A 2 u A , (A10) 

where A is the typical scale of perpendicular displacements of the field lines, and also the 
correlation length for such displacements. In our model, for typical modes, k±A ~ 1. It 
follows that k"j_i] tur b is equal to other estimates of the nonlinear interaction time scale. 

The maximum k± can be obtained from the condition that the energy transfer rate 
(see eq. (|A8|)) be greater than the resistive dissipation rate k\rj. We get 

kj < k ± , c l « ^. ( ™ . (All) 

Analogously, viscosity will be negligible as long as t~1 is greater than the viscous dissipation 
rate k\u±, which implies that 

^<t(ff> < A12 > 

where we have written the viscous diffusivity as v± to emphasize that it is the viscosity 
perpendicular to magnetic field lines which enters into this estimate. For a pure plasma the 
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viscous diffusivity may be large when the density and magnetic field are both low, but it 
will decrease as the field strength increases. 

In partially ionized gas, collisions with neutrals provide an additional source of 
viscosity. In this case the ambipolar diffusion coefficient r\ am b should be used instead of 77 in 
Eq. ([A For typical conditions in the cold interstellar medium (see Vishniac & Lazarian 



1998) 

Vamb w 5.3 x lO^n^N-jT^Bl cm 2 /s , (A13) 

where the letters with subscripts denote the normalized values of the parameters adopted. 
In other words, the electron density n is normalized on the density 0.03 cm -3 , neutral 
density N is normalized by 30 cm -3 , typical temperature T is assumed 100 K, while typical 
magnetic field 5 x 10 -6 Gauss. Such a high diffusivity implies k max l is only 10 3 for I = 30 pc, 
so that turbulence will be damped at ~ 0.03 pc scales in the cold interstellar medium. 

Rather high values of diffusivity follow from collisionless damping (see Foote & Kulsrud 
1969), but this only affects compressional modes in the cascade. GS95 showed that very 
little energy is transferred to the compressional modes from the shear modes on scales for 
which k± ^> Consequently, the suppression of compressional modes has little direct 
impact on the small scale structure of the magnetic field. 



The minimum k± covered by these scaling laws follows from equation (A"7). We find 

that 

VaY 



k±l > ■ (A14) 

Larger perpendicular scales will be described by weak Alfvenic turbulence. 

Are these scaling laws consistent with observations? The Sun provides our best example 
of strong MHD turbulence, but current observations of the solar corona and photosphere 
do not significantly constrain the power spectrum of the inertial range. On the other hand, 
the solar wind contains dynamical disturbances on a broad range of scales. Barnes (1979) 
has advocated that these should be interpreted as Alfven waves emitted from the solar 
corona and undergoing only weak dissipation in their outward journey. However Coleman 
(1968) has argued that the shorter wavelength disturbances are most naturally interpreted 
as a turbulent spectrum, and the solar wind shows evidence for local dissipation on small 
scales. A recent analysis by Leamon et al. (1998) supports the latter interpretation and 
argues that between 1.6 x 10~ 4 and 4 x 10~ 2 Hz perturbations in the solar wind are in the 
inertial range of strong MHD turbulence. (Quoted frequencies are measured at a spacecraft 
essentially stationary in a moving solar wind.) The power spectrum in this range has 
an index of roughly —0.7, consistent with a Kolmogorov index of —2/3. Since the solar 
wind advects disturbances over the spacecraft radially, while the unperturbed magnetic 
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field vector is strongly tilted by the Sun's rotation, this is essentially a measurement of 
the scaling law with respect to kj_. The observations also indicate that the turbulence is 
strongly anisotropic. This is in qualitative agreement with GS95, although a quantitative 
comparison would require an estimate of the fraction of the energy contained in nonlinear 
magnetosonic waves, and the local value of P ma gl ' P gas - This raises questions that have little 
impact on the nature of the small scale stochastic magnetic field structure, and we will not 
pursue them further. 

MHD turbulence in the ISM can be studied via scintillations. These measurements (see 
Spangler 1998) reveal that the density irregularities responsible for radio wave scintillation 
are elongated and aligned along the magnetic field. The scintillations themselves are 
believed to result from entropy fluctuations that arise from the mixing of plasma elements 
having different specific entropies. The entropy acts as a passive scalar and consequently 
its spectrum assumes the form of the turbulent energy spectrum (see Lesieur 1990). As the 
sharpest entropy gradients lie in the plane perpendicular to the magnetic field, we would 
expect to see the spectrum /c _5//3 . It is encouraging that the measurements (see Armstrong, 
Rickett & Spangler 1997) do get such a spectrum. Moreover, the elongation of fluctuations 
along the field lines follows from the disparity between h\ and k±. This may also serve as a 
hint in favor of Goldreich and Sridhar's modeQ 

Numerical simulations also seems to be consistent with the GS95 model. The high 
wavenumber resistive cutoff is given by Eq. ( [All ). Together with Eq. flA~7|) this sets the 
upper limit to k\\: 

We immediately see that, if b is the perturbation of magnetic field at the scale I, 

h 





( v y f \ 


b 


( v y /4 


— 

V A 


\IV A ) 




\IV A ) 



(A16) 



K'.L.c 

which is similar to the scaling relation^ reported in Matthaeus et al (1998). 



12 The range of scales over which the spectrum has this form is hotly debated. The emissivity of HI over 
the scale 10-100 pc has a shallower spectrum in the velocity space defined by galactic rotation (Lazarian 
1995). Lazarian & Pogosyan (1999) claim that this spectrum is consistent with the Goldreich & Sridhar's 
model of turbulence. However, here we are mostly interested in the smallest scales. 

13 Strictly speaking, in Matthaeus et al (1998) the relation k\\ iC /k oc b/B, where k — (k 2 c + k 2 ^ c ) 1 ^ 2 , was 
reported. The scaling above differs from this by ~ (b/B) 1 / 4 , which for the limited dynamical range available 
provides a difference of the order 30%. However, the difference may not be statistically significant. Peter 
Goldreich (private communication) attracted our attention to the fact that in Matthaeus et al (1998) fey 
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B. Previously explored ways to speed up reconnection 

B.l. Anomalous resistivity 

Anomalous resistivity is present in the reconnection layer when the field gradient 
is so sharp that the electron drift velocity is of the order of thermal velocity of ions 
u = [kT /m) 1 / 2 . For a detailed discussion of the physics see Parker (1979). The condition 
for appearance of anomalous resistivity is 

j > jcr = Neu . (Bl) 
If the current sheet has a width 5 with a change in the magnetic field AB then 

4*=^ . (B2) 

The effective resistivity increases nonlinearly as j becomes greater than j cr , thereby 
broadening the current sheet. We can therefore assume that j stays of the order of j cr , that 
is 

5*-^- . (B3) 

Expressing S through ion cyclotron radius r c = (muc) / \eB tot ) , where B tot is the total 
magnetic field (including any shared component) one gets 



. (V A \ 2 AB 

d Ri TV — — 



u ) Btot 

which agrees with Parker (1979) up to the factor AB/B tot , which is equal 1 in Parker's 
treatment. 

We can derive a reconnection speed by combining mass conservation and equation 
. We get 

t/ t/ rc ( Va \ 2 AB fan 

Vrec ~ V A — — — - . B5 

L x \u J B tot 



Equation (|E5|) shows that the enhanced reconnection velocity is much less than the 
Alfven velocity if the characteristic scale of the problem L x is much greater than the ion 
Larmor (cyclotron) radius. In general, "anomalous reconnection" is important when the 



is defined solely in terms of the background field and therefore if fcy is smaller than k±9, where 9 is the 
typical large scale bending angle, the local value of fen is not measured by this statistic. This will bias the 
measurements of the anisotropy towards a linear relationship. Definitely, new, more elaborate numerical 
tests are needed! 
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thickness of the reconnection layer in the Sweet-Parker reconnection scheme is less than S. 
This means that when 



S>(^L) ■ (BG) 



anomalous effects are important. The condition (|B6|) is easily fulfilled in the laboratory, 
where anomalous effects are widely observed. However, for typical interstellar magnetic 
fields the Larmor radius r c is ~ 10 7 cm. For interstellar Ohmic diffusivities rj w 10 9 cm 2 /sec 
and Va ~ u ~ 10 6 cm/sec we find that anomalous resistivity becomes important if 

S 2 V 

L x < « 10 11 cm « 0.01AU , (B7) 

V 

which is a very small scale for the ISM. All in all, anomalous effects are unlikely to be 
important for interstellar reconnection if the Sweet-Parker scheme is involved. 

In the presence of magnetic field turbulence the exact value of resistivity is much less 
significant. Nevertheless, anomalous resistivity may be important within the reconnection 
zone itself. By comparison with equation fl3T)|). we can see that anomalous effects will be 
important if 

or assuming Va ~ v\ and L x ~ / 



m 3/5 (Y±\ v * 



r <>'{ti) UfJ If) < B8 > 



/ < I0 13 cm « 1AU . (B9) 

Unfortunately, under normal conditions this is still a very small length scale. If I is less than 
~ 1 AU the speed of reconnection can be found from the condition V reCj i oca iX\\ ~ VaS, where 
Ay ~ k^ 1 is given by equation (|A7|) and we assume k± that enters this formulae ~ <5 _1 . 
Substituting 5 from equation (|B4|) we get 



Similarly, in the presence of the ambipolar diffusion in the gas with low degree of 
ionization the thickness of the reconnection layer becomes just barely greater than the ion 
Larmor radius (Vishniac & Lazarian 1998) so that we might expect the onset of anomalous 
effects, but even if they were marginally important they would not change the reconnection 
velocity significantly. In general we see that anomalous effects are of limited importance for 
reconnection in the interstellar medium. 

Recently this topic was revisited by Biskamp et al. (1997) and Shay et al. (1998) who 
showed that in a collisionless plasma the minimum width for the ion flow is actually the ion 
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skin depth c/u pi rather than the ion cyclotron radius. Since the former is just V^/ttj times 
the latter, this does not affect the conclusion that anomalous effects cannot explain fast 
reconnection in astrophysical plasmas 



B.2. Bohm diffusion 

Bohm diffusion is phrase which describes a process which is ubiquitous in laboratory 
plasmas, but which lacks a good theoretical explanation. Its characteristic feature is that 
ions appear to scatter about once per Larmor precession period. The resulting particle 
diffusion destroys the 'frozen-in' condition and allows significant larger magnetic field line 
diffusion. For an ion thermal velocity u ~ (kT/m) 1 ^ 2 and a cyclotron radius r c = muc/qB, 
where q is the ion charge, the Bohm particle diffusion coefficient is rjB ~ ur c « kT(c/qB). 
If we are concerned instead with the diffusion of field lines, we need to substitute V\ for 
kT/m so that r] Bohm ~ Va?V Since Va is of order u in the ISM this has no effect on our 
estimates. In the presence of Bohm diffusion, this coefficient should be used in place of 
Ohmic diffusivity (Parker 1979). However, we note that even if we make this substitution, 
it can produce fast reconnection, of order Va, only if r c ~ L x . It therefore fails as an 
explanation for fast reconnection for the same reason that anomalous resistivity does. The 
difference between the two processes is that anomalous resistivity vanishes when the scales 



involved become larger than the threshold values given by equations (|B7|) and ( p9|) , while 
we cannot define such a threshold for Bohm diffusion. 

We may repeat all our calculations in the presence of Bohm diffusion by substituting 
^Bohm instead of r] and as ryBohm ^> Vj ^ ne lower limit for the reconnection rates will increase 
by a substantial factor (^'Bohm/''?) 3 ^ 16 3> 1- 

A major shortcoming of this idea is that it is not at all clear that Bohm diffusion is as 
omnipresent under astrophysical conditions as in laboratory plasmas. Moreover, since our 
present study suggests that the actual reconnection rate does not depend on the value of 
diffusivity, this substitution is probably irrelevant. 



14 Shay et al. also found that the reconnection speed in their simulations was independent of L x , which 
would suggest that something like Petschek reconnection emerges in the collisionless regime. However, their 
dynamic range was small and the ion ejection velocity increased with L x , with maximum speeds approaching 
Va for their largest values of L x . Assuming that Va is an upper limit on ion ejection speeds we may expect 
a qualitative change in the scaling behavior of their simulations at slightly larger values of L x . 
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C. Tearing modes 

Tearing modes are a robust instability connected to the appearance of narrow current 
sheets QFurth, Killeen, fc Kosenbluth 1963| ). The resulting turbulence will broaden the 
reconnection layer and enhance the reconnection speed. Here we give an estimate of 
this effect and show that while it represents a significant enhancement of Sweet-Parker 
reconnection of laminar fields, it changes our results by only a small amount. One difficulty 
with many earlier studies of reconnection in the presence the tearing modes stemmed 
from the idealized two dimensional geometry assumed for reconnection ( [Furth, Killeen, 



fe Rosenbluth 1963|). In two dimensions tearing modes evolve via a stagnating non-linear 



stage related to the formation of magnetic islands. This leads to a turbulent reconnection 
zone ( |Matthaeus fc Lamkin 1985| ), but the current sheet remains narrow and its effects on 



the overall reconnection speed are unclear. This nonlinear stagnation stage does not emerge 
when realistic three dimensional configurations are considered. Indeed, it is easy to see 
that instead of islands one finds nonlinear Alfven waves in three dimensional reconnection 
layers. Therefore the tearing instability proceeds at high rates determined by its linear 
growth while the resulting magnetic structures propagate out of the reconnection region at 
the Alfven speed. 

The dominant mode will be the longest wavelength mode, whose growth rate will be 

The transverse spreading of the plasma in the reconnection layer will start to stabilize^ this 
mode when its growth rate is comparable to the transverse shear V4/A11 ( [Bulanov, Sakai,"%; 
Byrovatskii 1979| ). At this point we have 

Vreclocal « 7 A ~ V A ~. (C2) 

For a laminar magnetic field we have L x — An and these two equations imply 

rj \ 3/io 



v ™=Hvk) ■ (C3) 

This is substantially faster than the usual Sweet-Parker formula, but still quite slow. We 
note in passing that the value of A/Ay we assume here is below the minimum value proposed 
by Furth et al. (1963). However, it is consistent with the more accurate calculations of Van 
Hoven and Cross (1971). 



15 As a consequence of stabilization, tearing mode reconnection is stationary in a time-averaged sense. 
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If we consider instead a stochastic field with Aii <C L x , then we can express the local 
reconnection rate in terms of kn ~ AiT 1 . Substituting kn from equation ( |A7| ) and assuming 



k i ~ A we get 



while 



K 



reclocal 



and 



9/16 



kJ ~l.fJ It J ' (C6) 



The tearing instability broadens the reconnection region. For a completely ionized gas 
< k±^ max , and equation ( |C4D holds, provided that 

VAl >(^Y , (CT) 



which is nearly always true. In a partially ionized gas the tearing modes will be partly 
suppressed, and this mechanism for increasing V reCt i oca i will be relatively ineffective. 



Equation ( |G4| ) is the speed of reconnection for each flux element in a fully ionized 
gas, using Goldreich and Sridhar's model for the MHD turbulent cascade, and assuming 
the Sweet-Parker topology for the reconnection sheet. Equations (0) and (|C4] ) are not 
dramatically different. The exponent of the Lundquist number, TZi = (VaI/tj) changes by 
only 1/16. We also note that the frequency of reconnection, in this case the tearing mode 
growth rate, is ~ Va/X\\, so once more we can treat the field line structure as essentially 
static. 



D. Turbulence with an arbitrary spectrum 

Although in the main body of the paper we presented arguments in favor of Goldreich 
& Sridhar's (1995) model of turbulence, it is worthwhile considering reconnection when the 
relation between kn and k± is of the form 

h * r\kjy (^~y , (di) 

where p and m are unspecified positive constants. For GS p = 2/3, and m = 4/3. On the 
other hand, Matthaeus et al. (1998) have criticized the model of Goldreich and Sridhar, 
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partly on the basis of numerical simulations, and proposed an alternative model in which 
p = m = lQ. Using equation ( pi| ) instead of eq. (||), we get 



1(i-p)/p \v A 

which implies an outflow layer width of 

Lr\ 1/p ( vi \ m /P 



when I > L x and 



W-ltf) (£) , (D3) 

tf)* ~ (LJ,f» , (D4) 

when L x > I (compare to (|7|) and (||)). Therefore the upper limit on V rec imposed by large 
scale field line diffusion is given not by (0) but 



V rP , < V A min 



l 



vi) • (D5) 



which shows that for a generous range of p and m the outflow of fluid from the reconnection 
zone does significantly limit reconnection speeds. 

In order to recover limits on p or m which are consistent with fast reconnection we 
need to consider effects which might provide more stringent constraints on the reconnection 
speed. Some of our discussion cannot be reproduced for the general case without a model 
for the transfer of energy, as well as a scaling of eddy anisotropy. However, one important 
limit comes from considering the limits on small scale reconnection events. For the general 
case we replace equation (|13"D with 

Following the line of reasoning used in the main text requires that the probability a given 
magnetic field line will return back to the reconnection zone is small. Formally this means 
that p < 1. However, for p — 1 the local reconnection speed has already reached its maximal 
value Va, since in the case of locally isotropic eddies the reconnection geometry is basically 
the one proposed by Petschek (1964). For larger values of p reconnection proceeds at Va- 



16 Note, however, that the scaling of eddy anisotropy seen in their simulations is in agreement with GS95, 
once allowance is made for the existence of weak turbulence on large scales (cf. the discussion at the end of 
Appendix A of this paper). 
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Since p > 1 corresponds to a spectrum in which eddies become increasingly elongated across 
magnetic field lines on small scales, this choice is almost certainly unphysical as well. 

For p < 1 the Ohmic diffusivity provides the following estimate for the reconnection 
rate (compare with equation fll6|) 



V rec , global - k ]{ L x V rec , local -Va{—) - [ V J , (D7) 

which shows that the reconnection velocity does not depend on Lundquist number provided 
that p > 1/2. Otherwise Ohmic diffusion will limit the reconnection rate to the value given 
by eq. (|D7j). This will reach the Sweet-Parker reconnection rate when p = 0, that is, when 
the parallel wavelength does not increase as we go down the turbulent cascade. Physically, 
this is equivalent to assuming that the magnetic field lines are infinitely stiff, even when the 
local kinetic energy density is comparable to the magnetic field energy density. 

We note that even taking the relatively modest value of p = 1/2 implies that the 
magnetic field suppresses bending motions far more efficiently than one would suppose from 
simple energetic arguments. 

Finally, it is useful to see how these arguments would change if the Kraichnan spectrum 
turns out to be correct. We have earlier explained why we do not expect this to be the 
case, but this is the default model for many researchers. Also, it turns out that this model 
cannot be found as a special case of equation ( |D1|) . Instead, while most of the power is 
contained in modes with k± ~ k», the amplitude of the waves is given by 

^ « k- 1 ^- « k~\ (D8) 
u A V A 

which sets this case apart from the strongly nonlinear cascade, where the wave amplitude is 
~ kj 1 . If we combine the assumption of isotropy in wavevector space with equations (AT) 
and ( |A5| ) we find the usual result, for the Kraichnan spectrum, that 

v k ~Vi{kl)- x l*. (D9) 

Combining this with equation ( |D8|) we see that the characteristic displacement on a scale 
k^ 1 is 

< y2 > l/2_ ^_ I / m0 } 



This implies a field line diffusion, analogous to equation (pj) of 

dx \ P I \V A ) 
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which in turn leads to an outflow layer width, for L x < I, of 



(y 2 ) 1/2 ~i(¥) 2 (^T- (D12) 



/ ) \v A 

The upper limit on V re c,giobai set by field line diffusion is therefore 

WM 4 



V rec , global < V A -Z . (D13) 

To see that this is the actual speed of reconnection we need to evaluate V reCt i oca i- The 
Kraichnan spectrum will be truncated due to resistivity at a wavenumber given by 

Using the usual Sweet-Parker formula for the local reconnection speed we get 



2/3 / n \ 1/6 

V r ec,local ~ Va 



We can obtain the same estimate by ignoring the spectral cutoff, and calculating the 
appropriate value of k by requiring that ?7(|/ 2 )~ 1//2 equal the local Sweet-Parker reconnection 
speed. The upper limit on the global reconnection speed is recovered by multiplying V rec j oca i 
by kL x . We obtain 



1/2 

^rec, global ^ 



Comparing equations ( |D13| ) and ( |D16| ) we see that in most cases the reconnection speed 
will be given by equation ( P13j ). In other words, as usual, global field line diffusion provides 
the real constraint on magnetic reconnection in this case. 
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Fig. 1. — The geometry of magnetic field lines in 3D reconnection. The reconnected lines 
stretch and carry the conducting plasma with them. The plasma is also redistributed along 
the field lines. 
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Fig. 2. — a) The structure of the reconnection region when the field is turbulent. Local 
reconnection events happen on the small scale Ay rather than L x and this accelerates 
reconnection. The plasma is redistributed along the field lines in a layer of thickness (y 2 ) 1 ^ 2 , 
which is much thicker than the region ~ Aj_ from which the ejection of the magnetic field 
takes place, b) Local structure of magnetic field lines. 



